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Abstract 

In this study, the additional heat capacity which appear during 
the water dissociation of the proteins that are one of the soft materi- 
als, have been considered by the statistical mechanical methods. For 
this purpose, taking the electric field E and total dipole moment M as 
the thermodynamical variables and starting with the first law of ther- 
modynamics an equation which reveals the thermodynamical relation 
between the additional heat capacity in effective electric field ACe 
and the additional heat capacity at the constant total dipole moment 
ACm, has been obtained. It is found that, the difference between the 
heat capacities depends linearly on the temperature. To bring up the 
hydration effect during the folding and unfolding of the proteins the 
physical properties of the apolar dissociation have been used. In the 
model used for this purpose; the folding and unfolding of the proteins 
in the formed electric field medium have been established on this basis. 
In this study with the purpose of revealing the additional effect to the 
heat capacity, the partition functions for the proteins which have been 
calculated in single protein molecule approach by A. Bakk, J.S. Hoye 
and A. Hansen; Physica A, 304, (2002), 355-361 have been taken in 
order to obtain the free energy. In this way, the additional free energy 
has been related to the heat capacities. By calculating the heat ca- 
pacity in the effective electric field ACe theoretically and taking the 
heat capacity at constant total dipole moment ACm from the exper- 
imental data, the outcomes of the performed calculations have been 
investigated for Myoglobin and other proteins. 
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1 Introduction 



Now a days, a progress in the direction of investigation of the soft materials 
by the methods of statistical physics is realized. It is not possible to think a 
way of life without the soft materials; all biological structures, the molecules 
of the genetic code, proteins and membranes contain soft materials. Soft 
materials are the fundamentals of the life and an important component of 
the future technological developments. 

Proteins is a common name for the complex macromolecules which are 
formed by unification of a great number of amino-acids and which play 
unabandonable roles in the life times of all the living creatures. When 
considering the physical problem of formation of the structure of the natural 
proteins, one must keep in mind that proteins are quite different molecular 
objects. As a physical object the protein molecule is rather a big polymer 
composed of thousands of atoms and is such a soft material that from the 
physical point of view it is a macroscopic system. 

Every single atom forming the proteins occupies a place at a certain 
position as in the crystals. However, unlike to the crystals the position of 
each atom is unique to itself with respect to the neighboring atoms. Con- 
sequently, proteins represent a macroscopic system which is different from 
the crystals such that here the order is not periodical. A protein could be 
synthesized as a linear heterogeneous polymer which is formed by a genet- 
ically determined order where the remainders of different amino-acids are 
bounded to each other by peptide bounds. This polypeptide chain folds 
to a self unique conformation which is completely determined by the infor- 
mation in the ordering of the remainders of the amino-acids. The folding 
of the polypeptide chain to a three dimensional structure, principally is a 
reversible process and this depends on environmental conditions. In other 
words, thermodynamically, this process is passing from an unfolded state to 
another rather ordered macroscopic state. This case is named as the natural 
state of the protein molecule. 

2 Model 

For different aspects of the folding and unfolding transitions of the proteins 
various models have been proposed p] |2j . Simple but yet a notable model 
among then is the zipper model which defines the helix chain. In another 
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model, the idea of representing the solvent as dipoles in the unfolding of 
proteins was given by the work of Warshel and Levitt and later on in the 
applications of Rushell, Fan and Avbelj [31 HJ. In order to model the effect of 
the folding of the proteins Bakk asked for the inclusion of important physical 
properties in apolar dissolution. For this purpose he started by choosing the 
water molecules as classical dipoles (Sj. 



2.1 First Law of the Thermodynamics for the Proteins 

By taking the electric field E and the total dipole moment M as thermody- 
namical variables the first law of thermodynamics could be written for the 
proteins: 

dU = TdS - MdE. (1) 

Here, the electric field is not an external field but is used to model the ice- 
like behavior of the water molecules around non-polar surfaces. The electric 
field is a result of the effective behavior of the non-polar solvent applied in 
the unfolding of the protein. By considering the proteins as soft materials, 
the thermodynamical relation between the heat capacity at effective electric 
field AC E and the heat capacity at constant total dipole moment ACm has 
been put forward. Re-arranging Eq.(l) one could write 

dQ = dU + MdE. (2) 

To constitute such a relation, first of all the internal energy U is written 
as the total differential of the variables T and E and then substituted in 
Eq.(2) giving 

dQ = (^) E dT+[M + (^) T ]dE. (3) 

When the relevant heat capacities for constant M and E are written 
down in Eq. (3) one obtains: 

r)T7 f)E 

AC M = AC E + [M + (gi) T ](—) M (4) 

with the aim of expressing the term in the right hand side of Eq.(4) in 
terms of experimentally measurable quantities, Eq.(4) has been written in 
a different form, the total differentials of U(T,E) and S(T,E) have been 
calculated, the partial differentiations on both sides of equation have been 
performed and as a result the following expression has been found 

f) A/f B F 1 

AC M = AC E + T{—) E {—) M . (5) 
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Now, using properties of the total differential and also Maxwell relations 
constituted for the proteins 

AC M - AC E = -T{ 8 -^) T [(^) M f (6) 

has been obtained. 

Eq.(6) reveals the thermodynamical relation between the heat capacity 
at effective electric field ACg and the heat capacity at constant total dipole 
moment ACm- It is seen in the equation that; the difference of the heat 
capacities depend linearly on the temperature. ACe in Eq.(6) could easily 
be calculated theoretically and if ACm is measured experimentally, one 
could get a means of investigating the system in terms of thermodynamical 
quantities. 



2.2 Single Protein Approach 

Let us suppose that each atom has a dipole moment ~j? and interatomic 
interactions could be neglected. In the absence of an electric field each of 
the molecules has a continuous and randomly oriented dipole moment. 

On the non-denatured protein surface, in the positive bound charges 
(ions) of the water molecules an ice-like construction by the solvent effect 
becomes a matter of question. This in turn leads to the orientation of the 
dipole moments of the polar water molecules in the direction of field as if 
an external field were applied on the water molecules. This rigidness is for 
providing the water molecules to stay around the apolar surfaces. Thus this 
is not a real external field but reflects the characteristics of the solvation of 
the water molecules around the positive ions the non-polar surface in the 
protein folding. 

In this case the force on each dipole solely originates from the produced 
electric field. Since each atom experiences thermal agitation, the energy 
distributions of these atoms could be investigated by Classical Maxwell- 
Boltzmann Statistics [B]. The number of atoms dN in the energy internal 
U and U + dU is written as; 

dN = Cexp(-Uf3)dU (7) 

where (3 = 1/kT, k is the Boltzmann constant, T is the absolute temperature 
and C is a constant. U being the energy of the proteins, the energy gained 
by the dipole due to the applied field is given by 

U = -^-~e!. (8) 
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Let the angle between the axis of the dipole and the electric field be 9, 
then change in the amount of the electrical potential energy that the dipole 
possessed is obtained as 

dU =pEsm9d9. (9) 

When dU is substituted in Eq.(7), the number of atoms in the energy interval 
U and U + dU could be obtained in the following form 

dU = Cexp(p£ (3 cos 9)pEsin9d9. (10) 

If total number of atoms per unit volume is defined as N = Jq dN, then one 
writes for the total dipole moments as 

M = N < m >= / mcos9dN (11) 
J o 

where < m >= m which is the average dipole moment per atom in the 
direction of the formed electric field. When the equation expressing the 
total number of atoms is substituted in Eq.(ll) the solution in terms of the 
Langevin function is found as 

M = mL(a) (12) 

where pE(3 = a. 



2.3 Hydration Effect to the Single Protein Approach 

In addition to the energy coming from the external field, in order to introduce 
the pair interactions to the model pairing term must be added. For the 
purpose of obtaining an effective field in U e if two equations written initially 
are combined for a water molecule the effective field for this case becomes 

HE] 

U e = U DPE + U M f. (13) 

Substituting the mean-field solution of the hydration effect in the single 
protein approximation, one obtains for the dipole moment [1J Cj 

m= [ 27T d(f> r d9 cos 9 sin9 exp(-/3f7 e ). (14) 
Jo Jo 

Supposing that, distribution of the electric dipole moments has canonical 
distribution forms and making use of the fact that the orientations of p 
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vectors stay in the intervals 9 + d9 and 4> + d<J), from the definition of the 
solid angle d£l = sin 9d9d(ft, one gets the partition function for a protein as 

z = J J exp(-pU e )dn. (15) 

When solid angle dQ and the explicit expression for U e are substituted in 
the integral, one gets 

0bm 2 | 

z(T, E, N = 1) = 4vr exp — ~ — sinh (3e e . (16) 

By making use of the definition for the average dipole moment < m >= m/z; 
it is obtained in the form < m >= L(f3e e ). Here the effective energy is 
given bye e = e + 6<m>. The partition function of the system is the 
multiplication of N number of single particle partition functions. 

In thermodynamical investigation of the protein system, in order to in- 
troduce the chemical potential fi (here takes the place of pH), one has to 
take the system in a macro-canonical ensemble. In this study, the macro- 
canonical ensemble is obtained from the canonical ensemble. Here in addi- 
tion to be heat bath, it is supposed that the system is in a particle reservoir 
and if this reservoir is large, in the equilibrium, the particles reach an average 
value. 

In the first approximation, if the expression found in Eq.(16) for the 
canonical ensemble is extended for the macro-canonical case, without putting 
any limit on the particle number |Hj one could use the relation 



Z{T,E,v) = J2 ^[eM^HT,E,N = l)] N . (17) 

N=0 

The partition function of the undistinguishable particles which is given 
by Eq.(17) is obtained as 

Z(T, E, fi) = exp[exp(^)z(T, E, N = 1)] (18) 

where the expression £ = exp(/3/i) corresponds to the pH of the system and 
is the weight factor. In this case, 

is written for the energy of the system. Here, the expression $(T, E,/j>) = 
_\n_z(T^Ejj) - g ^ e thermodynamical potential. After the partition function 
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obtained by the substitution of Eq.(16) in Eq.(18) the thermodynamical 
quantities of the system are found as follows: 
Entropy is given by the formula 

S{T,E,») = -{^) E , IJl (20) 

which then leads to 



a(rp „ s , n r ,4exp^/x-^) 
S(T,E,fi) = A;[ln[exp( )] 

exp/3(/x - ^)[-47re e /3cosh(e e /3) 

+ We + 

27r(2 + /36m 2 )sinh(e e /3)] 



total dipole moment is given by the formula 



M(T,E,ri = -(—) T , lt (21) 



<)E" 

which then leads to 



x 47rexp(/9(/i-^-)) 
M(T,E,n) = P ^ 2 

(e e /3 cosh(e e /3) - sinh(e e /3)) 
and lastly the particle number is given by the formula 



N(T,E,h) = -(^) t ,e (22) 



leading to 

N(T,E,fi) 



4exp(/3(/x - ^-))^sinh(e e /3) 



3 Variation of the Additional Heat Capacity with 
Temperature 

In this study the variation of ACm with temperature has been investigated. 
In the macro-canonical ensemble the additional heat capacity in the effective 
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Figure 1: Graph of the variation of the additional heat capacities with tem- 
perature. 



electric field for the protein system is calculated from the expression 

ac -( du \ 1 ( dN \ u du ) i 2 m\ 

On the other hand, the variation of the internal energy with temperature is 
given by 

,dU. ,dU. .8U. ,8N. , nA . 

For (§?) T>E in Eq.(23) 

is substituted which is attained from Eq.(24). 

ACe is substituted in Eq.(6). M, S and iV that are necessary in Eq.(6) 
and Eq.(25) have been taken from Eq.(21), Eq.(20) and Eq.(22) respectively. 
The explicit forms of the above expressions are lengthy and thus are not 
written down. 

In order to compare ACm which is found theoretically with the exper- 
imental values of ACm they are shown on the same plot in Fig(l). The 
dashed curves represent the experimental study and the solid curves repre- 
sent the theoretical result. In this semi-phenomenological theory, the phys- 
ical quantities e , b and fi curves have been determined by fitting to the 
results of Privalov and Bakk . 

Red (Myoglobin), blue (Lysozyme), green (Cytochrome) and Yellow (Ri- 
bonuclease) solid curves represent the variation of the additional heat ca- 
pacity difference with temperature obtained from theoretical result. 



Table 1: Additional heat capacity values corresponding to transition tem- 
perature for four different proteins. 



Color 


Protein 


AHCD(^) 


Trans. Temp. (if) 


Red 


Myoglobin 


14.3 


301 


Blue 


Lysozyme 


9.3 


298 


Green 


Cytochrome 


8.6 


302 


Yellow 


Ribonuclease 


6.1 


320 



A H C D : Additional Heat Capacity Differences 



4 Conclusions 

In this study a theory of the heat capacity depending on electrostatic field 
and constant dipole moment in the transition between the fold and unfold 
states of the protein molecule according to Bakk J3| model, has been pre- 
sented. The foresight of the performed calculations have been discussed by 
taking into account the experimental studies of Privalov and Makhatadze 
on Myoglobin, Lysozyme, Cytochrome and Ribonuclease in the interval 
275 - 350K [SI El- 

Biological and other macromolecular systems exhibit great variations in 
heat capacity with respect to temperature. The unfolding of the globular 
proteins in water is a classical example. Change in the environmental con- 
ditions such as temperature, increase or decrease in pH, interaction of the 
solution dissolvers with protein molecule groups, gives rise to the change in 
the structure of the proteins. 

When the conditions return to their initial state these molecules return 
to their former structures. The folding and then unfolding of these molecules 
occur by a thermodynamical mechanism. When the graph of the variation 
of the additional heat capacities with temperature is examined one could 
state that the behavior of the proteins exhibits approximately an universal 
behavior. The apex in the heat capacity curve indicates the denaturation 
transition. For each of the four proteins; additional heat capacity difference 
and transition temperatures Tt correspond to different values and the nu- 
merical values have been presented in Table 1. For every protein in native 
state the heat capacity is less than the value in the denature field. The 
transition between the native and denature state is a first order transition. 

An additional reason of the observed increase in the heat capacity is the 
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hydration effect. The transformation of the apolar solvents into water, takes 
place with a considerable increase in the heat capacity. A conclusion drawn 
here is that; depending on the unfolding of the proteins, the additional 
increase in the heat capacity is a consequence of the contact of the internal 
groups during the unfolding of the proteins. The hydration effect originates 
from the low solubility of the apolar compounds in water. 

One of the characteristics of the folding and unfolding of the proteins 
is that, the temperature interval of the process depends on environmental 
conditions and particularly on the pH values of the solution. This task is 
supplied by the term ( = exp(/?/x) in the partition function. Variations in 
the enthalpy difference relevant to pH will be presented in a forthcoming 
article. 
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